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ABSTRACT 

We present an infrared study of a z = 0.872 cluster, SpARCS J161314+564930, with the primary aim 
of distinguishing the dynamical histories of spectroscopically confirmed star-forming members to assess the 
role of cluster environment. We utilize deep MIPS imaging and a mass-limited sample of 85 spectroscopic 
members to identify 16 24 /im-bright sources within the cluster, and measure their 24 /im star formation rates 
(SFRs) down to ~ 6 M© yr _1 . Based on their line-of-sight velocities and stellar ages, MIPS cluster members 
appear to be an infalling population that was recently accreted from the field with minimal environmental 
dependency on their star formation. However, we identify a double-sequenced distribution of star-forming 
galaxies amongst the members, with one branch exhibiting declining specific SFRs with mass. The members 
along this sub-main sequence contain spectral features suggestive of passive galaxies. Using caustic diagrams, 
we kinematically identify these galaxies as a virialized and/or backsplash population. Moreover, we find a mix 
of dynamical histories at all projected radii, indicating that standard definitions of environment (i.e., radius and 
density) are contaminated with recently accreted interlopers, which could contribute to a lack of environmental 
trends for star-forming galaxies. A cleaner narrative of their dynamical past begins to unfold when using a 
proxy for accretion histories through profiles of constant (r/r2oo) x (Av/cr^); galaxies accreted at earlier 
times possess lower values of (r/r 200) x (Av/a v ) with minimal contamination from the distinct infalling 
population. Therefore, adopting a time-averaged definition for density (as traced by accretion histories) rather 
than an instantaneous density yields a depressed specific SFR within the dynamical cluster core. 
Subject headings: Galaxies: clusters: general - Galaxies: clusters: individual: SpARCS J161314+564930 - 
Galaxies: high-redshift -Galaxies: star formation - Infrared: galaxies 



1. INTRODUCTION 

A hallmark of cos mology has been the es tablishment of hi- 
erarchical formation (IWhite & FrenklH~991l) in which growth 
of structure in the Universe proceeds via a "bottom-up" sce- 
nario: matter condenses into low mass haloes which eventu- 
ally merge to form larger structures, culminating with galaxy 
clusters. Clusters subsequently evolve through the accretion 
of galaxies and groups along cosmic filaments. This continual 
build-up of the cluster gives rise to distinct populations: a viri- 
alized component of older galaxies, and a younger population 
that was recently accreted from the surroundin g low-density 
field dBalogh et al.lll998k lEllingson et afll2001h . As clusters 
are thought to be hostile environments, member galaxies are 
exposed to various mechanisms that could potentially sup- 
press their star formation , including ram-press ure stripping 
(e.g., iGunn & Gotj 119721) , strangulation (e.g., lLarson et al.l 
1980), and ga laxy harassment, perhaps preceded by an ini - 
tial burst (e.g.. lBarnes & Hernquisll99lHMoore et al.ll 19981) . 
These mechanisms should induce marked differences in the 
young and in-situ population, as the latter has long endured 
the harsh conditions of the dense cluster. 

Indeed, the local environment in which a galaxy resides 
is known to strongly correlate with several galaxy proper- 
ties, such as star formation rate, stellar mass, color, and mor- 
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phology. Extensive observational efforts at low redshift have 
yielded a paradigm for galaxy dependencies, such that the 
densest reg ions at z ~ are devoid of star forma tion ac- 
tivity (e.g.. iGomez et alJl2003l: IBalogh et alJl2004al) and in- 
stead harbor massive (e.g., Kauffmann et al.ll2004|). red (e. 
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early-type passive galaxies (lDresslerlll980l) . 

Recently, the nature of these correlations at higher red- 
shift, z ~ 1, has become a rather contentious issue. While 
it is well established that higher redshift clusters contain in- 
creased star f ormation activity compared to their local coun- 
terparts (e.g.jButcher & Oemlerlll978l 119841 lEllingson etal] 
1200 U ISaintonge et alj 120081 Webb et al. in preparation) in 
parallel to the rapid decl ine of the cosmic star formation rate 
(SFR) since z ~ 1 (e g..lLillv et aljfl996l iMadau et all [l996l 
iLe Floc'h et al.ll2005l) . there has not yet been a clear conver- 
gence on the SFR-density relation at this epoch. Some studies 
have observed a reversal of the z ~ relation such that the 
predominant sites of star formation have migrated to denser 
regions by z ~ 1 (lElbaz et al.l 120071: ICooper et al.l 120081: 
iLi etalJ201ll) . while many other groups fin d the local relation 
is already in place in z ~ 1 clusters (e.g.. | Patel et al] 120091 
l2lfiltlKov7rna et alJUoTollMuzzin et alj|2012l Webb et al. in 
preparation). These discrepancies most likely stem from dif- 
fering selection effects (stellar mass versus luminosity limited 
samples) and the varying degrees of densities probed (cluster 
versus group environments). 

Moreover, the situation is muddled by the interdependence 
between environment and stellar mass — whether these corre- 
lations with environment are causal or incidental is still am- 
biguous, since massive galaxies preferentially reside in dense 
regions. Indeed, many properties seem to be also governed 
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by stellar mass, i ncludi ng color (e.g., iBell & de Jong! 1200 It 
Kauffmann et alj 12003b and SFR (e.g.. iBrinchmann et alj 
2004tlNoeske et al.ll2007l) . As surveys have attempted to un- 



tangle this covariance between stellar mass and environment, 
a host of trends have been revealed to depend unilaterally on 
one, or equally on each property. For instance, both mass and 
environment have a separable effect on the fraction of star 
forming galaxies as seen by a decline in star-forming galax- 
ies with increasing de nsity and mass while t he other parame- 
ter remains fixed (e . g..lPoggianti et al. | 2008t iPeng et aLll201 Ob 
ISobral et al.l 1201 It iMuzzin et al.1 12012IK In contrast, many 
studies have found that the specific star formation rate (SSFR; 
the star formation rate per stellar mass) for star-forming g alax- 
ies is correlated with stellar mass in different environments 
but f ails to exhibit any dependence on the local environment 
(e.g. iKauffmann et alJ2004tlPoggianti et alJ2008tlPeng et alJ 
I2010t ILu et al.M2012t IMuzzin et alJl2012h . A possible expla- 
nation for this discordance that has been adopted by many 
groups is a rapid quenching timescale that suddenly alters 
the host galaxy's star-forming classificatio n (e.g., color) pre- 
ceding a decline in the observed SS FR dPeng et al.l 120101: 
IMuzzin et alJIMTllWeTzeTetiniMll) . 

Alternatively, a flat trend in SSFR with environment could 
suffer from radial projection effects. If radial/density bins har- 
bor a mixture of both physically high and low radius galax- 
ies, any correlation with star forming properties could poten- 
tially get washed out in projection space. Moreover, this could 
preferentially affect star forming galaxies over the quiescent 
population if they are inherently more biased towards radial 
co ntamination. Thi s is similar to the rationale put forward 
bv lPatel et al.l (1201 ll) . who attributed a declining SFR-density 
(over all quiescent and star-forming galaxies) at 0.6 < z < 
0.9 to a combination of two effects: a varying fraction of pas- 
sive and active galaxies with density, and suppressed SFRs at 
higher densities. 

The recent work of lHaines et al.l (120121) provides further in- 
sight into this contamination scenario via the accretion histo- 
ries of cluster galaxies. Utilizing caustic dia grams from the 
Millennium Simulations, lHaines et alj d2012|) effectively iso- 
late active galactic nuclei in projected velocity/radius space 
and determine that they are primarily an infalling population 
in spite of their low projected radii. Applying the results of 
accretion histories from their work could potentially incite 
the emergence of an alternative environmental trend for star- 
forming galaxies that relies on more of a time-averaged rather 
than an instantaneous definition of density. 

Here, we present an infrared study of a z = 0.872 
galaxy cluster drawn from the Spitzer Adaptation of the 
Red-sequence Cluster Survey (SpARCS: IWilson et al.l 120091: 
IMuzzin et aT] 120091: IDemarco etail 120101). We utilize ex- 
tensive optical spectroscopy ( IMuzzin et all 1201 2|) and deep 
Spitzer-MlPS 24 /im observations to pinpoint dusty star- 
forming cluster galaxies. By adapting a kinematic approach 
to the star formation histories of these galaxies via caustic 
diagrams, we intend to alleviate some of the confusion in en- 
vironmental trends caused by radial projection effects. 

The paper is outlined as follows. In if2] we present the 
SpARCS survey and 24 /im observations. We briefly describe 
our computation of stellar masses and star formation rates 
in $3] Our results are presented in SjU including stellar age, 
kinematic, and star formation trends. We discuss the implica- 
tions of these trends in |5] an d introduce a kinematic approach 
to classifying the star formation histories of cluster galaxies, 
along with a new interpretation of the environmental depen- 



dence on the SSFR of star-forming galaxies. We conclude 
in Throughout the paper we assume a cosmology with 
H = 7010ns- 1 Mpc~\ tl M = 0.3, = 0.7. Stellar 
masses and SFR s are based on a Chabrier initial mass function 
(IChabrierl2003l) . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. The SpARCS Sur\>ey 

The Spitzer Adaptation of the Red-sequence Cluster Sur- 
vey (SpARCS) is a ~ 45 sq. deg. survey with deep z'-band 
imaging from the CFHT and CTIO, designed to produce 
a large, homogeneously selected sample of massive galaxy 
clusters at z > 1. By combining the z'-passband obser- 
vations with IRAC imaging from the Spitzer Wide-area In- 
fraRed Extragalactic (SWIRE) survey, SpARCS has discov- 
ered high-z massi ve cluster candidates using either the red- 
sequence method dGladders & Yedl2000|: IMuzzin et a .1120081) 
or the stellar bump sequence method (IMuzzin et al.l 120131) . 
With ~ 200 massive cluster candidates, including ~ 12 z > 
1 spectroscopically confirmed clusters dWilson et all 120091 
IMuzzin et alj|2009t IDemarco et al.ll2010l) , SpARCS currently 
has one of the largest repositories of z > 1 galaxy clusters. 

2.2. SpARCS J161314+564930and the GCLASS Cluster 
Sample 

This work presents an infrared study of SpARCS 
J1613 14+564930, a rich galaxy cluster at z = 0.872 
from the G emini Cluster Astro physics Spectroscopic Survey 
( GCL AS S : IMuzzin et al]|20 1 21) that was discovered using the 
red-sequence method. It is an extremely massive cluster with: 
a velocity dispersion of 1350 ± 100 kms -1 ; r2oo, the radius 
at which the cluster density is 200 times the critical density, 
of 2.1 ± 0.2Mpc; and M 2 oo = 26.1±|1 x 10 14 M (Wilson 
et al. in preparation). This cluster field has ~180 spectro- 
scopic redshifts from GMOS-N on Gemini. Galaxies within 
3375 km s -1 of the cluster velocity dispersion (2.5(7,,) are 
considered to belong to the cluster, yielding 95 confirmed 
cluster members. As GCLASS is a 3.6 /un-selected sur- 
vey, the spectroscopic data is rest-frame i/-band selected and 
therefore close to a stellar mass-selected sample. The final 
sample includes 85 cluster galaxi es ab ove our mass complete- 
ness limit of 2.0 x 10 9 M Q (see ffOt . 



2.3. Spitzer-MIPS Imaging 

Our primary data set derives from the Multiba nd Imaging 
Photometer for Spitzer ( MIPS: iRieke et alj|2004l) aboard the 
Spitzer Space Telescope ( I Werner et al.l 120041) at 24 fj,m with 
the aim to detect any dusty emission associated with the spec- 
troscopically confirmed members. With a 5.4 sq. arcmin field 
of view, MIPS provides sufficient coverage of the entire spec- 
troscopic area. The observations were part of the Guaran- 
teed Time Observer program and completed in 2008 (proposal 
ID 50161) with an exposure of 1200 seconds per pixel. The 
MIPS image was reduced using a combination of the Spitzer 
Science Center's MOPEX software and a suite of IDL rou- 
tines we developed to further optimize background subtrac- 
tion (Webb et al. in preparation will provide more details). 

2.4. Source Detection and Photometry 

A source catalog of the entire MI PS field contai ns flux den- 
sities estimated with DAOPHOT dStetsonlll987h . and posi- 
tions using a photometry pipeline ( lYedll991l) . We plot the 
24 fim differential number counts for our catalog in Figure [TJ 
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they are in good agreement with the published 24 fim counts 
from iPapovich et all (120041) . We determine the catalog com- 
pleteness limits by locating the flux at which our counts devi - 
ate by > 2a from the expected value o flPapovich et al.1 (120041) . 
which occurs at fluxes below ^70 /iJy. This corresponds to an 
infrared luminosity of — 3 x 10 10 L Q at z — 0.872. 



S 



10' 



C/3 



* 1 1 



? 



X Papovich et al. (2004) 
■ This work 



100 



1000 



S v (uJy) 



FIG. 1 . — The 24 /im differential number counts of our catalog, normalized 
to the Euclidean slope (filled squares). We compare to published counts from 
Papovich et al. 1 2004) to determine a completeness limit of 70 /Jy. 



3. ANALYSIS 

3.1. Counterpart Identification 

We search for optical counterparts within ^2" to locate 
cluster members with 24 fim emission; we detect 16 MIPS 
cluster members which constitute the focus of our study. To 
quantify the number of chance alignments between 24 /im and 
optical sources, we drop 2" apertures at random positions on 
the MIPS field and count the number of times we find one 
or more source. We perform 100 realizations at each flux 
level within 40-400 /iJy in 10 fiiy steps. At 70 fiJy, the com- 
pleteness depth of our MIPS catalog, we expect 5% of MIPS 
counterparts to be falsely identified. This translates into one 
chance alignment in our counterpart sample of 16 MIPS clus- 
ter members. 

In order to assess the effect of the cluster environment on its 
constituent galaxies, it is essential to have a control sample of 
field galaxies that are distinct from the cluster. Our spectro- 
scopic campaign for the entire GCLASS survey successfully 
obtained redshifts for a significant number of foreground and 
background sources in the fields of the clusters with identi- 
cal selection effects as c luster galaxies; thes e sources consti- 
tute our field sample (see lMuzzin et al.ll2012l) . We employ the 
same criteria to search for MIPS counterpart emission to the 
field and cluster galaxies, but restrict the field redshift range 
to 0.84 < z < 1.0 in order to probe the same epoch as the 
z = 0.872 cluster. This leaves us with ten 24 /im-bright field 
galaxies above our flux limit. 

3.2. Stellar Masses 

We compute stellar mas ses using a technique s imilar to 
iBell & de Jongl d2001l) and iKauffmann et alj ( 120031) . In par- 
ticular, we utilize a spectral feature that provides informa- 
tion on the age of the stellar population but is fairly insen- 
sitive to dust, the 4000-A break. Using the lBruzual & Charlotl 



( 120031) stellar population synthesis models with solar metal- 
licity, assuming a Chabrier IMF, and adopting a star forma- 
tion timescale of r = 0.3 Gyr, we infer a M/L ratio at 3.6 fim 
from the strength of the 4000-A break. We than convert to a 
stellar mas s using the measured 3.6 fim luminosity (for more 
details see iMuzzin et al.l 120121) . In order to obtain a mass- 
limited sample, we do not include sources below our mass 
completeness limit of 2 x 10 9 M© (IMuzzin et al.ll2012l) . 

3.3. 24 fim Star Formation Rates 

The mid-infrared luminosity of a galaxy probes thermal 
emission from dust grains that has been reprocessed from UV 
light, and therefore offers a clean measurement of the star 
formatio n rate (SFR) tha t is minimally affected by dust ex- 
tinction (lKennicuttl[T9 98). Moreover, rest-frame 24 fim traces 
dust heated by younger stars an d therefore provides a mea- 
sure of the instantaneous SFR (ICalzetti et al.ll200"7l) . How- 
ever, at z = 0.872, 24 fim corresponds to 13 fim rest-frame, 
which has recently been shown to contain a higher contribu- 
tion of dust heated from intermediate age stars, which is ev- 
idence for star for mation over longer timescales (~l-2Gyrs; 
Sali m et al.ll2009l) . As there i s a rapid decline i n the cosmic 
star formation since z ~ 1 (Lilly et al.1 1 19961 iMadau eTaT] 
Il996t LeFloc'het al.1 120051) . this could result in an overes- 
timate of the sta r formation as derived fr om observed 24 fim 
flux. In parallel, Rodighi ero et al.l (120101) found that extrapo- 
lating an infrared luminosity from the monochromatic 24 fim 
flux leads to an underestimate by a factor of ~ 1.6 compared 
to Herschel measurements observed at 100 and 160 fim for 
0.5 < z < 1.0 galaxies. Despite some potential systematic 
effects in the 24 /im-derived SFRs, the results from this study 
are minimally affected, as we rely so lely on relative d iffer- 
ences within our own sample (see also lPatel et al.ll20TTl) . 

We calculate the 24 /im-derived SFR by converting the 
MIPS flux into a total infrared l uminosity using an average 
of the ICharv & Elbazl (120011) and iDale & Helovl(l2002l) tem- 
plates. We employ the relation in iKennicuttld 1998b to calcu- 
late a SFR and convert to a Chabrier-IMF based SFR using a 
factor of 1.65. Based on the depth of our MIPS images, we 
probe down to ~ 6 M Q yr _1 . 

3.4. Active Galactic Nucleus Contamination 

We assume that the MIPS sources are dominated by star 
formation, with little or no contamination from active galactic 
nuclei (AGN); this seems reasonable given that only ~ 5% of 
infrared field galaxies at z ~ 1 have their total IR luminos- 
ity dominated by AGN emission dFadda et al.ll2010l) . and the 
contribution declines drastically with low mid-infrared flux. 
In fact, AGNs only begin to dominate the infrared emission 
for 24 fim flux es above 1 .2 mJ y at z ~ 0.8, and are missing 
below 0.8 mJy dFu et al.ll2010l) . We note that our sample con- 
tains only two sources at the high flux end, at 0.83 mJy and 
0.91 mJy, both of which are still below t he level where AGNs 
prevail. Moreover. [Martini et al.l (120091) measured the occur- 
rence of AGNs in cluster galaxies and found that while the 
AGN fraction increases with redshift, it is still only 1.5% at 
z = 0.81. 

Regardless, we have performed two quick checks for ob- 
vious AGN sources. We first compare the equivalent width 
of t he [Oil] line and the d epth of the 4000-A break, follow- 
ing IStasihska et alj ( 120061) . Our sample contains ten MIPS 
sources with measurable [Oil] emission, all of which are 
all consistent with the normal star forming galaxies shown 



4 



Noble et al. 



in fig. 7 from IStasinska et al.l (120061). We also exploit in- 
frared color diagnostics (e.g., lLacv et ail 120041: ISaiina et"aT] 
2001) using data from the Infrared Array Camera (IRAC; 
Fazio et all2004h . All but two MIPS cluster members display 
infrared colors consistent with stellar- and PAH-dominated 
sources (i.e., star-forming) at z ~ 0.9. The other two sources 
lie on the edge of the region encompassed by AGN, but 
we note there is a high contamination from PAH-dominated 
sources at this location; these are the same two sources de- 
scribed previously with high mid-infrared fluxes. In the ab- 
sence of any concrete evidence for pure AGNs, we assume 
the MIPS flux is dominated by star formation. 

4. RESULTS 

This analysis presents an infrared study of a single 
GCLASS cluster at z = 0.872. As such, we focus pri- 
marily on the properties of the MIPS population and utilize 
the extensive spectroscopy for verification of MIPS cluster 
members hip and an es t imate of stellar age. We refer the 
reader to iMuzzin et al.l (120121) for a detailed description of 
the spectroscopic selection criteria and completeness, as well 
as a comprehensive analysis of the GCLASS clusters as seen 
through optical spectroscopic measures, for example, [Oil] 
star-formation rates. 

4.1. Age as a Function of Stellar Mass 

Here we investigate the depth of the 4000-A break (D4000) 
as a function of stellar mass for the MlPS-detected clus- 
ter members compared to the larger sample of spectroscopic 
members. The 4000-A break, defined as the ratio of the inte- 
grated flux density at 4 000 A-4100A to th at blue-ward of the 
break, 3850 A-3950 A dBalogh et alj 1999ft . arises from an ac- 
cumulation of absorption lines and increases in depth for old 
and metal-rich stellar populations. It can therefore be used as 
a proxy for stellar age. In Figure [2] we plot the depth of the 
4000-A break for all cluster galaxies (filled circles), as well 
as a field population over the same epoch, 0.84 < z < 1.0 
(open squares). We highlight cluster and field galaxies de- 
tected at 24 fim with green diamonds. The MIPS cluster 
galaxies are primarily coincident with the coeval field pop- 
ulation, spanning the same range of stellar mass and D4000. 
A Kolmogorov-Smirnov test between the MIPS cluster mem- 
bers and general field population reveals a 76% chance they 
derive from the same parent population, in stark contrast to 
that of cluster members without MIPS detections, where the 
null hypothesis is rejected at a high significance of 0.02%. 
Compared to the spectroscopically-confirmed members not 
detected at 24 (im, MIPS members form a young envelope 
of the cluster population at a given mass, which suggests they 
were recently accreted from the field. 

We also distinguish between galaxies with [Oil] emission 
and those without (i.e., below ~ 1 (3) A equivalent width for 
the highest (lowest) signal-to-noise spectra); in Figure|2]cyan 
circles have measurable [Oil] which is indicative of ongoing 
star formation, while purple circles represent quiescent (or 
extremely dust-enshrouded) galaxies. The majority of MIPS 
galaxies are [Oil] -detected in the cluster and field, albeit at a 
slightly reduced fraction in the cluster (70 and 90%, respec- 
tively). This is not surprising given the recent conclusion from 
Webb et al. (in preparation) that optical studies are not signifi- 
cantly biased against dusty star formation. We define a D4000 
limit of 1 .4 (dotted line), above which > 90% of cluster mem- 
bers without [Oil] emission reside, signifying a minimum age 



for the older, quiescent cluster galaxies. We note that a vari- 
able D4000 as a function of stellar mass does not significantly 
affect our conclusions, and therefore adopt the simpler flat de- 
lineation. In all subsequent plots, we differentiate galaxies 
by their nominal age as traced by the strength of the 4000-A 
break, rather than [Oil] detections as in Figure [2] 
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FIG. 2. — The strength of the 4000-A break as a function of stellar mass 
for the cluster (z = 0.872; filled circles) and coeval field (0.84 < z < 1.0; 
open squares) population. The galaxies with [Oil] emission (star forming) 
are plotted in cyan, while those without any [Oil] emission are purple. A 
green diamond denotes coincident 24 fim emission; there are three MIPS de- 
tections, all with [Oil] emission, at the location of the thickest green diamond 
with a D4000 of 1.3 and stellar mass of ~ 1.7 X 10 9 M Q . The dashed ver- 
tical line illustrates our mass limit of 2.0 X 10 9 Mq, and the dotted hori- 
zontal line indicates our cut in D4000 in subsequent plots. We note that the 
two most massive cluster members have both been identified as the brightest 
cluster galaxies, although only one is detected at 24 fim. 



4.2. Relative Velocity as a Function of Radius 

In the left panel of Figure [3] we plot the relative velocity 
of all cluster members as a function of clustercentric radius, 
as defined by the projected distance to the brightest cluster 
galaxy (BCG), the brightest 3 .6 fim cluster memb er, which 
is also the most massive (see iLidman et alJ 120121 for a de- 
tailed discussion of the BCG sele ction in SpARCS, includ- 
ing this cluster). As discussed in i l4.11 we now separate the 
galaxies into old (red circles) and young (blue circles) clus- 
ter populations based on D4000, and highlight 24 /im-bright 
members with green diamonds. Barring two MIPS galaxies 
at low projected velocities and radii (both of which are clas- 
sified as older), MIPS sources tend to display larger veloci- 
ties on average (|Av| = 1422 ± 205.8kms _1 ), in contrast 
to a tighter distribution expected from a virialized population 
(\Av\ = 995.1 ± 89.60kms _1 for older galaxies without a 
MIPS detection). Therefore, the velocity-radial space encom- 
passed by MIPS galaxies further supports the idea that they 
have been recently accreted. 

The histogram in the right panel of Figure[3]emphasizes this 
point: the older cluster population (red hashed histogram) dis- 
plays a Gaussian distribution of velocities and primarily falls 
within the cluster velocity dispersion (1350kms _1 ), while 
the young (blue hashed histogram) and MIPS (solid green his- 
togram) galaxy populations exhibit bimodal and/or flat dis- 
tributions, peaking at velocities greater than ±1000kms _1 . 
Moreover, the younger members within the cluster avoid the 
central velocity bin completely. 
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Given that 38% of MIPS cluster galaxies have a steep 4000- 
A break with D4000 > 1.4 (red circles with green diamonds), 
we might expect to see some overlap in their velocities com- 
pared to those in the larger spectroscopic sample of older 
galaxies (i.e., the confirmed members that are not detected 
by MIPS; red circles without a green diamond). However, a 
Kolmogorov-Smirnov test between the velocities of all MIPS 
galaxies (green diamonds) and older cluster galaxies (the 
black histogram shows old cluster members with old MIPS 
sources removed) rejects the null hypothesis with marginal 
significance — there is only a 6% probability that the two dis- 
tributions derive from the same parent population. It therefore 
seems likely that the MIPS galaxies trace a younger, recently 
accreted population that is not yet virialized with the cluster. 
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FIG. 3. — The relative velocity as a function of clustercentric radius for 
all spectroscopic members. Red circles corresponds to older galaxies with 
D4000 > 1.4 and blue circles denote younger galaxies, D4000 < 1.4. The 
histograms correspond to the number of galaxies within 900 km s~ 1 velocity 
bins; the colors correspond to the same types of galaxies in the legend on the 
left. The black histogram represents the population of old galaxies with the 
MIPS sources removed (6 in total). The green histogram shows all MIPS 
sources, regardless of age. 



4.3. Environmental Dependence on the Specific Star 
Formation Rate 

Given the implication that MIPS galaxies represent a 
younger, infalling population of cluster members, we might 
expect to detect a correlation between the specific star for- 
mation rate (SSFR; star formation rate per stellar mass) of 
MIPS galaxies and local environment. If the majority of MIPS 
galaxies do in fact belong to an infall population, there should 
be more star formation in the cluster outskirts compared to the 
cluster core, where galaxies have had their star formation shut 
off as they fall deep into the cluster. 

In Figure|4] we plot the SSFR as a function of local environ- 
ment, with binned averages shown in the lower panels. The 
upper panels show the SSFR for every MIPS member (green 
stars), as well as an upper limit for every spectroscopic mem- 
ber, determined from the 3cr completeness limit of the SSFR 
for a given stellar mass. Errors on the averages are determined 
from 100 bootstrap resamplings of the data in each bin. 

We use two different proxies for environment in order to 
compare to previous studies, and to confirm we are not sub- 
stantially biased with either parameterization. In the left pan- 
els we show clustercentric radius, where the core of the cluster 
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is defined by the location of the BCG The right panels con- 
tain local galaxy density calculated from the distance to the 
10th nearest neighbor. Although in the literature these two 
parameters are often used interchangeably, they sample differ- 
ent mechanisms: the former is a better indicator of the global 
environment and cluster potential, whereas the local density 
has a proximate effect on galaxies. However, both parameters 
still play a significant role in the properties of cluster galaxies, 
and we theref ore investigate trends with each separately (see 
iLi et al.ir2012l for a detailed discussion on the two methods). 
While both measurements suffer from projection effects, the 
benefit of using local density is that it is straight-forward to 
correct each galaxy fo r spectroscopic com pleteness given its 
stellar mass and radius dMuzzin et al.l2012h . which essentially 
provides us with a 100% complete sampling of density. We 
note that our radial measurement does not require a complete- 
ness correction as th e completeness bias for mass is similar at 
all radii (see fig. 4 in lMuzzin et al.l20T2l) . To compute density, 
each galaxy is first given a weight based on its completeness. 
We then interpolate the distance to the 10th nearest neighbor 
for each source by summing the weighted values for galaxies 
above our mass limit of 2.0 x 10 9 M Q , and compute the den- 
sity as £ = 10/ (ndf ). The drawback of this method is that 
we lose information on smaller density scales. 

Nevertheless, our two environment parameterizations con- 
vey the same trend: the average SSFR of star-forming galax- 
ies (i.e., MIPS galaxies; green stars) is mostly independent of 
environment (see lower panels of Figure H). A best-fit line to 
the average MIPS SSFR in radial and density bins reveals a 
slope consistent with zero at the level of 1.6er (—0.29 ± 0.18) 
and 1.2cr (—0.23 ± 0.20), respectively. In Figure|4] we show 
the best fit to the normalization of the SSFR with the slope 
fixed at zero. Though seemingly surprising, this flat trend is 
consistent with other cluster surveys over many redshifts and 
using various star formation trace rs, including UV studies at 
0.16 < z < 0.36 dLu et al.ll2012l) and optical star forma tion 
indicators at z ~ 1 dPeng et alJl20iotlMuzzin eTai1l2012l) . 

However, our MIPS sample is SFR-limited, and thus not 
complete in SSFR for any given mass except at the highest 
SSFR values, log(SSFR) > -9. The flat trend in SSFR might 
be the result of skimming the high SSFR galaxies off the top 
of a deeper correlation. For example, if we assume that more 
massive galaxies reside at lower radii, probing deeper in SFR 
could preferentially bring down the average SSFR at low radii 
compared to outer radii, uncovering a non-flat SSFR trend. 
We check for this effect by assigning a 3a upper limit for the 
SSFR of cluster members without a MIPS detection (black ar- 
rows in upper panels of Figure @J and include these limits in 
the average SSFR (black circles in lower panels). The upper 
panels of Figure|4]reveal that there is no underlying mass seg- 
regation in the radial/density properties in the population not 
detected by MIPS (i.e., the limits are flat with environment), 
which suggests that a lower SSFR-limit would not necessar- 
ily result in a non-flat trend for star-forming galaxies galaxies. 
Moreover, the slope of the average SSFR assuming all cluster 
members have some star formation flattens out even further 
with density (—0.08 ± 0.10, dotted black line in lower right 
panel). 

If we instead treat the galaxies without MIPS detections 
as quiescent (i.e., non-star forming) and investigate the inte- 
grated SFR (above our SFR-limit) per total unit stellar mass 
as a function of environment (cyan triangles), namely the total 
MIPS SFR divided by the total stellar mass of all cluster mem- 
bers in each bin, we find a decreasing SSFR with increasing 
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(decreasing) density (radius). We note that both spectroscopic 
and MIPS members have the same target selection, so they 
have the same completeness rates as a function of mass and 
radius. The integrated SFR per total stellar mass is a proxy 
for the number of star forming galaxies compared to quies- 
cent galaxies at each radius and reveals a lower fraction at 
denser environments. A fit to the integrated SSFR versus den- 
sity yields a slope of —0.61 ± 0.19 (dashed cyan line), and is 
therefore inconsistent with zero at the > 3ct level. 

This depression of star formati on in the cluster core is in ac- 
cordance with local clusters (e.g. JKauffmann et al.ll2004l) and 
suggests that the SSFR-density relation is already established 
in the highest-density regions at z <~ 0.9. This is in contrast 
to z ~ 1 field studies that found a reversal in th e SFR-density 
relation dElbaz et al.ll2007l:ICooper et al.ll2008l) , although they 
were limited to lower density environments. Instead, we com- 
pare our results to a simil ar IR study of a z = 0.834 clus- 
ter from iPatel et all (120091) who probe a wide-range of envi- 
ronments and utilize a 24 fim stacking analysis on all cluster 
members to determine SSFRs. They uncover a trend of de- 
creasing SSFR with increasing density (open squares in the 
right panel of Figure |4}, analogous to our total SSFR. This 
relation prevails even when controlling for stellar mass: the 
open, teal squares correspond to lower masses of 2.0 x 10 10 < 
M/M Q < 6.3 x 10 10 , and the open orange squares represent 
all galaxies with M > 6.3 x 10 10 M Q in the lPatel et alf(l2009l) 
sample. We note the normalization of our trend is higher, 
but we probe down to an order of magnitude lower in stellar 
mass (2.0 x 10 9 Mq); as lo wer-m ass systems typically display 
slightly higher SSFRs (see tH.4| below), a higher SSFR in our 
data is expected. 

In Figure|4]we also plot the SSFR required for stellar mass 
to double (assuming a constant SFR) by z = (dot-dashed 
horizontal line). We note that the integrated SSFR per total 
stellar mass (cyan triangles) is consistent with or below this 
limit in all environments, signifying that a majority of clus- 
ter members (including quiescent galaxies) have already ex- 
perienced most of their mass growth. However, the average 
MIPS SSFR (and average SSFR for all members with lim- 
its) lie above this line: the star-forming galaxies have yet to 
undergo the bulk of their activity, assuming they can sustain 
their SFR. Perhaps we are witnessing the primal growth stage 
of these galaxies, before they have properly assimilated into 
the cluster and been exposed to any environmental quench- 
ing. Alternatively, if any of these galaxies have already expe- 
rienced a suppression of their star formation, the quenching 
timescale must be rapid enough that any environmental trend 
with SSFR gets washed out: SSFRs are immediately termi- 
nated, precluding an intermediate phase. In other words, why 
is there a lack of variation in the star-forming population de- 
spite a changing fraction of star-forming galaxies with envi- 
ronment? We will return to this flat environmental trend in |J5] 
and provide an alternative explanation. 

4.4. The Correlation between Star Formation Rate and 
Stellar Mass 

We can investigate whether any MIPS members deviate 
from their expected SFRs and SSFRs given the tight corre- 
lation of increasing SFR with stellar mass for star-forming 
galaxies, which retains only 0.2 dex scatter at 2 ~ 1. Re- 
cent studies have observed this star-forming main sequence in 
the field out to z ~ 2, which monotonically shifts to higher 
SFRs with increasing red shift (jNoeske e t al.l2007llElbaz et alj 
120071: iDaddi et all 120071) . In Figure |5] we plot the main se- 



quence trend in the field at various redshifts: z = 0.1 from the 
Sloan Digit al Sky Survey (SDS S; lBrinchmann et al.ll2004b as 
analyzed bv lElbaz et alj d2007l) : z ~ 1 fro m the Great Obser - 
vatories Origins Deep Sur vey (GOODS: lElbaz et al.l 12007b : 
and z ~ 2 from GOODS (IDaddi et al.l 120071) : we have con- 
verted each trend to Chabrier-IMF based masses and SFRs for 
consistency with our own sample. We also plot our 24 /im- 
derived SFRs for MIPS cluster members at z ~ 0.9 along 
with our field sample over 0.84 < z < 1.0. Due to system- 
atic differences in SFRs, stellar masses, and selection criteria 
(e.g. mass-limited versus luminosity-limited samples), we re- 
frain from quantitative comparisons between our trend and the 
field samples from GOODS and SDSS. Our own field sample 
offers a more suitable comparison, and in fact, displays a sim- 
ilar trend as the larger GOODS sample at z ~ 1, albeit with 
slightly higher SFRs. 

Immediately obvious in Figure|5]is the distinction between 
the old and young MIPS populations: the majority of massive, 
older members (filled red circles) lie well below the expected 
correlation, while the younger members (filled blue circles) 
tend to follow the main sequence trend observed in the z ~ 1 
GOODS field (solid black line) and the GCLASS field (open 
squares). Moreover, this double-sequence in SFRs is unique 
to the cluster as there is only one field galaxy in our sample 
with a significantly low SFR for its stellar mass. 

The double-branched distribution is similarly manifested 
in the SSFR, plotted in the lower panel of Figure The 
least massive galaxies, which are inherently younger (see Fig- 
ure O, have the highest SSFRs, while more massive galaxies 
have substantially lower SSFRs. While t his is consistent with 
the evolution of the mass function (e.g., Kodama et alj|2004b 
which shows that the massive end of the galaxy mass function 
in clusters is in place by z ~ 1 and the evolution between 
2 = — 1 consists of a build-up of the < 10 11 M Q end, it 
could also be a completeness effect as we do not probe the 
region of low mass and low SSFR (Figure[5J. 

In hopes of identifying further differences, we split the 
24 /im members into two cases based on their proximity to 
the expected SFR for their given mass at z <~ 1, as shown by 
the dot-dashed pink line in Figure [5] This line corresponds 
to log(SSFR)= —9.3. We designate the population above the 
line as the main sequence population, and those below the line 
as the sub-main sequence given they are in fact closer to the 
z ~ field trend and lie in a substantially disparate space 
from the GCLASS z ~ 1 field sample. 

In Figure [6] we stack the spectra in each population sepa- 
rately (nine and seven members in the main sequence and sub- 
main sequence groups, respectively) with an inverse weight- 
ing based on the spectroscopic completeness. Specifically, we 
investigate the equivalent widths of the [Oil] doublet (3727 A) 
and a Balmer absorption line, H5 at 4100 A (see Table[T|); their 
relative strengths impart a timescale of star formation, with 
strong [Oil] indicating current activity and HS representing 
more prolonged star formation. The striking contrast in the 
stacked spectra hints at differences in the star formation his- 
tories of these two populations. The main sequence group 
has moderate to strong [Oil] emission, relatively deep HS ab- 
sorption, and strong Balmer absorption features. Based on 
the spectral classifications of Dre ssier et al.l (11999b . these are 
most likely dust-obscured galaxies undergoing bursty star for- 
mation, e(a) galaxies. They also border on the classification 
of e(c) galaxies, which experience normal, continuous star 
formation. As these main sequence galaxies are coeval with 
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FIG. 4. — Upper panels — The individual SSFRs for all cluster members as a function of clustercentric radius (left panel) and density using the 10th nearest 
neighbor (right panel). The green stars correspond to MIPS cluster members, and the black circles with arrows are 3cr upper limits on the SSFR given the 
stellar mass for the cluster members without MIPS detections. Lower panels — The SSFR in bins of clustercentric radius (left) and density (right). The green 
stars correspond to the average SSFR of detected MIPS members in each bin, i.e., the SFR divided solely by the mass of MIPS members and averaged over the 
number of MIPS galaxies in each bin. Black circles represent the average SSFR for all cluster members, assuming a 3cr upper limit on the SSFR for undetected 
MIPS members. Cyan triangles denote the integrated SFR (above our SFR-limit) per total unit stellar mass: a sum of the total 24um-derived SFR divided by 
the total stellar mass of all spectroscopic members in each radial bin, which essentially probes the fraction of star-forming galaxies. The horizontal solid green 
line illustrates the best fit line to the green stars with the slope fixed to zero, i.e., the normalization of the average MIPS SSFR. The dot-dashed horizontal line 
corresponds to the required SSFR for the mass to double by z = 0. The two remaini ng lines in the lowe r right panel depict the best linear fits to the corresponding 
binned values with density. The open squares represent the stacked IR SSFRs from Patel et al. (2009) in a z = 0.834 cluster, where teal squares correspond to 
galaxies with masses of 2.0 X 10 10 < M/Mq< 6.3 X 10 10 , and the grey squares represent all galaxies with M > 6.3 X 10 10 Mq. 



z ~ 1 GCLASS field, and follow the expected SFR versus 
stellar mass trend, they have probably been recently accreted 
from the field, and therefore belong to the infall population. 

On the other hand, the sub-main sequence population has 
weak [Oil] emission and a strong 4000-A break, marked by 
a sudden drop at [Call] K and H and lines. The shallow 
H<5 absorption is also indicative of older systems: as massive 
stars die, the H6 absorption fades. The lack of strong Balmer 
absorption precludes the possibility of poststarburst systems. 
The stacked spectrum for this population is more consistent 
with that of a passive galaxy with only a slight hint of star 
formation; it is classified as right on the border of an older 
k-type spectrum and e(c) galaxy, as there is some [Oil] emis- 
sion present. We list equivalent widths for [Oil] and H6, along 
with D4000 in Table Q] 

We also note the possibility that galaxies on the sub-main 
sequence could be AGNs, rather than simply star-forming 
galaxies with lower star formation activity. This would also 
be an interesting explanation as it could imply that the AGN 



TABLE 1 

Spectroscopic measurements from the weighted, stacked 

SPECTRA. 





Sub-main sequence 


Main sequence 


Oil EW 


5.4 


22.0 


US EW 


2.4 


4.7 


D4000 


1.46 


1.20 



is responsible for qu enching the star fo rmation through feed- 
back. For example, Page et al] ( 120121) have claimed to find 
evidence for quasar-mode feedback in the form of suppressed 
star for mation rates for the mo st luminous AGNs at z = 1 — 3 
(but see lHarrison et al.l(l2012l) for an alternative result). In this 
scenario, we could be witnessing the residual star formation 
in AGNs that are still actively accreting and just beginning to 
peak in luminosity. 

5. DISCUSSION 
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FIG. 5. — Top — SFR versus stellar mass for MIPS cluster galaxies, sep- 
arated in color by D4000. Open squares are the field population from 
0.84 < z < 1.0. The dotted, solid, and dashed black lines correspond to 
the field trends at 2 = 0, 1, 2, respectively (converted to a Chabrier IMF). 
The pink dot-dashed line denotes the boundary between the two populations 
of the cluster galaxies: main-sequence vs sub-main sequence, as defined by 
their star formation rate. The vertical and horizontal gray dashed lines indi- 
cate our mass completeness, M > 2.0 X 10 9 Mq, and SFR detection limit, 
respectively. Bottom — SSFR versus stellar mass. Symbols and lines are the 
same as the top panel. 

We have analyzed a single z = 0.872 massive cluster from 
the GCLASS sample from an infrared perspective, investigat- 
ing how dusty star-forming galaxies behave in the broader 
context of all spectroscopically-confirmed cluster members. 
We discuss our results here, and merge them with a detailed 
kinematic analysis of each MIPS member. This allows us to 
access the role of environment on star-forming galaxies utiliz- 
ing their previous exposure to high densities, rather than their 
instantaneous environment. 

The ages and stellar masses (Figure |2) of the MIPS mem- 
bers compared to both cluster and field galaxies suggests that 
they have been recently accreted from the field. The line- 
of-sight velocity distribution for MIPS members (Figure |3) 
exhibits analogous results, as they hug an outer envelope of 
velocities consistent with the younger cluster members and 
what would be expected from an infall population. MIPS clus- 
ter members have been observed to share certain properties 
with field galaxies in previous studies, but a re not simply a 
parallel population in all cases. For example. iKocevski et alj 



FIG. 6. — The resulting st acke d spectra of MIPS cluster members from 
each population identified in i|4.4l main sequence population (cyan spectrum) 
and sub-main sequence population (pink spectrum). Individual spectra are 
weighted by their spectroscopic completeness, as determined by their stellar 
mass and clustercentric radius. 



(1201 ll) determined that although MIPS sources have a spatial 
distribution within the cluster typical of an infalling popula- 
tion, their spectral properties reveal burstier episodes of star 
formation than their field counterparts at z ~ 0.9, perhaps 
provoked by harassment and mergers during assembly. While 
the stacked spectrum of our main-sequence galaxies is con- 
sistent with e(a) type galaxies — possibly indicative of bursty 
star formation — we do not see evidence for enhanced activity 
compared to the field. Moreover, we uncover an additional 
branch of 24 /im members that exhibit older, more quiescent 
spectral features, suggesting the existence of an environmen- 
tal quenching mechanism that could occur following a possi- 
ble (though not definite) initial triggering of star formation. 

We find the SSFR of MIPS members to be independent 
of the projected local cluster environment: they maintain the 
same level of star formation at all densities and radii (Fig- 
ure |4]i._Althou£h tlrisjsTOnsi^ (e.g., 
iPeng et alJl2()10l:lLu et afl2012tlMuzzin et alJl2012l) . it is sur- 
prising given an observed increase in the fraction of star- 
forming galaxies with decreasing density over various star- 
forma tion indicators at z < 1, for example, fO IIl emission 
lines dPoggianti et alj 120081: iMuzzin et all 120121) : Ha lumi- 
nosity dSobral et alj|201 lb : and 24 fim flux (e.g.. IPinn et alj 
1201(1 and the work presented here as shown by the cyan 
triangles in Figure |4j. Many studies have invoked rapid 
quenching timescales to explain the disparate trends, such 
that star-forming galaxies are never observed in an interme- 
diate stage — they are either active or completely quenched. 
However, when we investigate the SFRs and SSFRs as a func- 
tion of stellar mass (Figure |5J, we discover two distinct star- 
forming branches: one that is in line with the expected main 
sequence in the z <~ 1 field and undergoing obscured star 
formation, and another that displays depressed levels of star 
formation for a given mass and has a stacked spectrum con- 
sistent an older k-type galaxy. Moreover, these sub-main se- 
quence galaxies are on average more massive; it is the less 
massive galaxies that are forming the bulk of the stars (though 
this could be a selection effect). This is consistent with a 
downsizing trend in mass assembly such that star formation 
shifts to less massive galaxies as the Un i verse ages (see als o 
iKodama et aT] [20041 iFeuiner et all 120071 ISobral et alJUoTlh . 
Perhaps this sub-main sequence is a population of MIPS 
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galaxies that have reached higher density regions or even the 
cluster core sometime in their past and have been exposed to 
some sort of environmental quenching; however, this seems 
inconsistent with the fact that no environmental dependence 
on SSFR is observed (Figure |4}. Either quenching is ubiqui- 
tous throughout the cluster with little dependence on the en- 
vironment or we are not measuring the environment properly. 
This raises the question: are we truly sampling MIPS mem- 
bers in the cluster core? We address this issue in the subse- 
quent sections. 

5.1. Velocity Distributions of Two Star-Forming Populations 

In Figure [7] we re-analyze the line-of-sight velocities as a 
function of projected radius for the MIPS members, account- 
ing for the double-sequence in SFRs. We highlight the MIPS 
sources on the main sequence with cyan squares, and those 
populating the sub-main sequence with pink squares. There 
seems to be a rather clear distinction between the two pop- 
ulations, with the majority of sub-main sequence galaxies 
possessing lower velocities, while the main sequence galax- 
ies avoid velocities below 1000 km s -1 ; with a few excep- 
tions, the main sequence MIPS members adhere to an outer 
envelope in radial-velocity space. More importantly, at any 
given projected radius, we observe a mix of main sequence 
and sub-main sequence star-forming galaxies. If we return 
to the notion that the SSFR has no significant dependence on 
environment, this result offers a hint of an alternative inter- 
pretation: it indicates that low radial bins are contaminated 
by high-velocity galaxies, which could be infalling galaxies 
(outside the cluster core) that have fallen far enough into the 
gravi tational potential w ell to attain large line-of-sight veloc- 
ities dHaines et al.ll20l2b . Some of these systems may in fact 
be at physically low radii, however, they could be on their 
initial (or even a subsequent) pass through the cluster and 
therefore are still distinct from a virialized core population 
that has spent a much greater amount of time in dense envi- 
ronments. As these systems have retained a nominal level of 
star formation for their stellar mass, they augment the average 
SSFR in low radial bins, or in what is referred to as the cluster 
core. The antithesis contamination occurs at large projected 
radii where the average SSFR is pulled down from low veloc- 
ity interlopers — possibly galaxies at large physical radii that 
have already felt the effects of the cluster environment and 
therefore are not representative of an infalling population. In 
tandem, these two effects yield an overall flat SSFR with pro- 
jected radii and densities. 

We attempt to identify dyna mically distinct popula tions in 
Figure [7] utilizing the results of iMahajan et al.1 (1201 ll) . where 
they use the dark m atter hydrodynamical simulations from 
iBorgani et al.l (120041) to extract a sample of 93 mock clusters 
from simulations of 117 haloes, 105 of which have M200 > 
1.4 x 1O 14 M0. From these clusters, they statistically quan- 
tify the fraction of virialized, backsplash, and infall galax- 
ies at various projected radii and velocities. These classifi- 
cations correspond to decreasing timescales of accretion, re- 
spectively: virialized galaxies were accreted when the clus- 
ter core was forming and are now passively evolving; back- 
splash galaxies have had time to pass through the cluster core 
and experience the first effects of the cluster environment, but 
have not yet amalgamated with the virialized population and 
therefore represent an intermediate phase of accretion (e.g., 
iBalogh et alj l2000h ; infall galaxies were recently accreted 
from the field and have not reached pericenter in their orbit 
around the cluster. We plot the delineations in velocity-radius 



phase space in Figure |7j the corresponding fractions of each 
classification, along with the applications to our entire spec- 
troscopic sample are summarized in the right-hand panel. 

We assign each galaxy in our spectroscopic sample as viri- 
alized, infalling, or backsplash based on its age and, for the 
case of MlPS-detected members, its SFR. For cluster mem- 
bers that are not detected by MIPS, we simply assume that 
older members are virialized and younger galaxies have been 
recently accreted (i.e., infalling). For MIPS galaxies, we as- 
sume any galaxy on the main sequence belongs to the infalling 
population, while those on the sub-main sequence are virial- 
ized if they are older, or are part of the backsplash population 
if they are younger. 

The consistency b etween the simulated fractions from 
iMahaian et al.l d201 lb and our cluster sample, using this sim- 
plistic criteria, is quite good. For example, the solid outlined 
box in Figure [TJrepresents galaxies within 0.5r2oo and l.Oert,, 
where 89%, 8%, and 3% of galaxies are expected to belong 
to the virialized, infall, and backsplash populations, respec- 
tively; we have 52 cluster members within this area, therefore 
the predicted fractions correspond to 46.3 virialized, 4.2 in- 
falling, and 1.6 backsplash members in our sample (the sec- 
ond column in the right panel displays the expected num- 
ber within our sample in each box). Of the seven MIPS- 
detected galaxies within this region, three fall on the star- 
forming main sequence (cyan squares), which suggests they 
are recently accreted members that have retained their gas 
and star formation — this number is just below the expected 
number of infalling galaxies from the simulations (4.2), and 
increases to seven after we include the members not detected 
by MIPS. The remaining MIPS members in this region lie on 
the sub-main sequence (pink squares), and therefore display 
lower SFRs. Additionally, there is a clear divide in the age 
of these four galaxies, with two containing shallow 4000 A 
breaks (filled blue circles). The depressed SFRs of these two 
galaxies, along with their younger stellar populations, suggest 
they are backsplash galaxies which were stripped of their gas 
reservoir as they passed through the cluster core, but have not 
had enough time to completely virialize; this is also consistent 
with the expected backsplash fraction (two galaxies in total). 
The last two MIPS members (pink squares surrounding red 
circles) are both older and have suppressed SFRs, suggestive 
of passive, virialized galaxies which statistically comprise the 
bulk of this region. Including the galaxies without MIPS de- 
tections that have older ages, our classification scheme yields 
43 virialized galaxies in this region, which is consistent with 
the number from the simulations (46.3). We note there is 
(seemingly) one MIPS outlier at both large projected radius 
and velocity that has a low SFR and an older stellar popula- 
tion. However, this region is still expected to contain 40% 
virialized galaxies, so it is not unfeasible for less active sys- 
tems to inhabit this phase space. 

We check the robustness of our classifications with respect 
to the depth of the MIPS data by assuming the SFR-limit is 
a factor of two lower (3 M Q yr^ 1 ), and all galaxies currently 
not detected by MIPS have this SFR. We then simplistically 
assign them to the above classifications based on their stel- 
lar mass and which star-forming sequence they populate. We 
find that two infalling galaxies would switch to a backsplash 
classification and three virialized galaxies would be infalling; 
therefore, there is no significant change in the numbers as- 
suming a deeper SFR-limit. 
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FIG. 7. — Left — The absolute line-of-sight velocity as a function of projected radius. MIPS members belonging to the main (sub-main) sequence of SFR are 
highlighted with cyan (pink) squa res. The horizontal lines correspond to cr v and 2a v , and the vertical line represents 0.5r200- The boxed areas correspond to the 
divisions in Mahaian et al. 1 201 if) and are identical in both panels. Right — The first number on the left in each rectangle correspond s to the fraction of gal axies 
classified as virialized (orange), infall (purple), and backsplash (gray) in projected phase space bins from dark matter simulations in Mahaian et al. l2uTl|). We 
then apply the fractions to the total number of spectroscopically confirmed cluster members in each bin and list it after the colon. The third number represents 
the actual number of each type of galaxy within that phase space based on its age and SFR (for MlPS-detected members). We assume a simple classification in 
which: virialized galaxies consist of non-star forming older galaxies (red circles) and MlPS-detected older galaxies on the sub-main sequence (red circles with 
pink squares); infalling galaxies include non-star forming younger galaxies (blue circles) and MlPS-detected galaxies on the main sequence, both young (blue 
circles with cyan squares) and old (red circles with cyan squares); and backsplash galaxies comprise younger MlPS-detected galaxies on the sub-main sequence 
(blue circles with pink squares). The last column in each bin (green number in parentheses) represents the number of MIPS sources of that galaxy type that are 
included in the actual number counts based on the classification scheme stated above. 



5.2. Utilizing Caustic Diagrams as Accretion History 
Predictors 

The mixing of different dynamical histories at all projected 
radii casts a doubt on the interpretation that environment has 
no effect on the SSFR of star-forming galaxies. However, our 
attempt to classify each galaxy as infalling or belonging to the 
cluster core is circumstantial. 

A complementary method to kinematically differentiate be- 
tween galaxies arises with the use of caustic profiles (i.e., 
constant velocity-radial lines) in physical radial phase space, 
where inf alling and vir ialize d cluster galaxie s become dis- 
tinct (e.g.jMamon et all|2004t iGill et all2005l:lMahajan et all 
1201 U lHaines et alj|2012l) . Moreover, this is inherently linked 
to the epoch of accretion, which illustrates the importance of 
a dynamically-defined environment: the time-averaged den- 
sity a galaxy has been subjected to is more significant than its 
immediate s urroundings. 

Recently, lHaines et alJ (120 12h . stacked 30 clusters in the 
Millennium Simulation at z = 0.21 and investigated the rela- 
tion between projected me asurables and the ac cretion history 
of the cluster (see fig. 3 in lHaines et alJ|20l2T) . Galaxies that 
were accreted when the cluster was first forming have a large 
dynamic range of velocities at the smallest projected radii, but 
only occupy a narrow silver of velocities (Av/v a ~ 0) over 
larger projected radii. Conversely, recently accreted galax- 
ies permeate all projected velocities and radii, but are pri- 
marily concentrated along and outside a caustic line. These 
trumpet-shaped caustic profiles roughly correspond to lines 
of constant r x v in projected phase space. Recently accreted 
galaxies are therefore the dominant contaminant compared to 
galaxies that were accreted earlier, specifically in low radial 
bins. Moreover, star forming studies are biased towards pick- 
ing out these contaminants, assuming recently accreted galax- 
ies have managed to retain their activity and are classified as 
star-forming. This could explain why SSFR trends utilizing 



the entire galaxy population (quiescent and star-forming) or 
stacked SSFRs p reserve a correlation with environment (e.g., 
iPatel et aT1l201 lb . while those solely for star-forming galaxies 
flatten out. 

Motivated by the possibility of identifying distinct accre- 
tion histories within our own sample via caustic profile di- 
agnostics, in Figure [8] we plot the strength of the 4000-A 
break as a function of (r/f2oo) x {Av/a v ), since constant 
values of this parameter correspond to caustic profiles which 
trace out accretion epochs. There is a clear movement to- 
wards higher values of (r/r 200) X (Av/a v ) for shallower 
breaks, meaning lower values on the x-axis correspond to 
older galaxies that were accreted at earlier times. We calcu- 
late a linear Pearson correlation coefficient of —0.41 between 
the parameters, with a high significance of > 99.99%. To 
ensure this trend is not driven purely by the correlation be- 
tween D4000 and stellar mass, we have separated the cluster 
members into two mass bins, above and below 3 x 10 10 M Q , 
which is roughly the median mass of the sample. The relation 
persists for both mass bins, meaning that for cluster galax- 
ies of the same stellar mass, there is a progression towards 
larger (r/r 200) x (Av/cr v ) values with decreasing time since 
the last burst of star formation. Assuming higher values of 
(r/r 200) x (Av/a v ) trace galaxies that were more recently 
accreted (which is analogous to fig. 3 in lHaines et al.ll2012l) . 
the trend in Figure [8]indicates the cluster environment has an 
effect on its constituent galaxies. 

In order to isolate contamination in low projected ra- 
dial bins from infalling/backsplash interlopers, we divide the 
space into three distinct regions based on the location of main 
and sub-main sequ ence galaxies. Combining the results of 
lHaines et al.l (120121) and [Mahaian et alJ ( 1201 11) . these regions 
should primarily correspond to: virialized galaxies that were 
accreted when the cluster core was forming; a mix of all types 
but where backsplash galaxies are most likely to exist; and 
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infalling galaxies that were recently accreted — at low, inter- 
mediate, and high values of (r/r2oo) x (Av/a v ), respec- 
tively. We calculate the median D4000 in each region for both 
mass bins (filled stars), which confirms the declining age (i.e., 
epoch of accretion) with (V/r 2 oo) X (Av/a v ). 
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FIG. 8. — The strength of the 4000-A break versus (r/r 2 oo) x (An/(T„) 
for all cluster members, divided in two mass bins. The MIPS members on the 
main and sub-main sequence are highlighted with cyan and pink squares, re- 
spectively. The solid gray line shows the best fit line for all galaxies. The me- 
dian D4000 values for both mass bins in each region (denoted by the vertical 
dashed lines) are plotted as filled stars and indicate a declining accretionary 
sequence with (r / r2oo ) x ( /cr v )• The error bars are taken from the stan- 
dard deviation on the bootstrapped median in each bin. We have labeled the 
regions based on accretion times. We note the three galaxies with the lowest 
caustic values are artificially plotted at (r/r2oo) x (At)/ff») = 0.006 to 
reduce the size of the plot, but in fact have even lower values. 

We demonstrate the utility of these caustic regions in ref- 
erence to our own sample in Figure [9] We plot lines of con- 
stant r x v at 0.1 and 0.4 (r2oo x &v), based on the regions 
in Figure [8] The innermost caustic line seems to preferen- 
tially select the older (virialized) population of cluster mem- 
bers that were most likely accreted at early times, including 
three out of the seven MIPS sources belonging to the sub- 
main sequence branch. Although rare, it is not inconceivable 
that recent s tar formation has occ urred within the virialized 
population; (Mah aian et al.l l20l"ll) estimate that 11% of the 
virialized population are galaxies with ongoing or recent effi- 
cient star formation that could be attributable to rapid flybys. 
However, we note that their level of star formation is certainly 
suppressed compared to the main sequence. 

Between the two caustic profiles, we expect a distribu- 
tion of galaxies that were accreted early, recently, and some- 
where in-between. Indeed, in this intermediate region of 
0.1 < r/r-200 x v/a v < 0.4, there exists a mix of old, main 
sequence galaxies and young, sub-main sequence galaxies, 
which could result from a backsplash population. The sig- 
nificance of the backsplash po pulation has been de monstrated 
before in various simulations. Bal ogh et al.l (12000) found that 
over half of the galaxies within l-2r2oo have been inside the 
virial radius at an earlier time and rebounded outward. More- 
over, 90% of these backsplash galaxies have plunged deep 
into t he potential we ll, within the inner 50% of the virial ra- 
dius dGill et al.1120051) . 

The region exterior to the 0.4 (r2oo x cr„) caustic profile 
should preferen tially pick out gal axies that have been re- 
cently accreted (lHaines et al.l I2.Q12D . Based on the fractions 



from iMahaian et al.l (1201 ll) . it should addi tionally favor in- 
falling over backsplash galaxies. Moreover. iGill et ail (120051) 
observed that basksplash galaxies are kinematically distinct 
from infalling galaxies, as the latter have much higher rel- 
ative velocities. Although this velocity bimodality can get 
slightly washed out in projected space, there still exists an av- 
erage trend toward lower velocities for backsplash galaxies. 
Notably, this region contains the majority of young main se- 
quence galaxies, which are most consistent with an in fall pop- 
ulation given their age and SFRs. As stated in c]5.1| the lone 
sub-main sequence galaxy in this region could statistically be 
a virialized galaxy outside the core; however, we also note 
that we expect one chance al ignm ent between MIPS sources 
and spectroscopic members ( 33.11 . which could be manifested 
here. Alternatively, this galaxy and another old, sub-main se- 
quence galaxy at high projected radius could both belong to 
groups that have been accreted into the cluster. In this sce- 
nario, the lower SFRs of these galaxies could be due to pre- 
processing (e.g., IZabludoff & Mulchaevlll998 ) in the group 
environment prior to accretion. iMcGee et al.l (|2009) found 
that the accretion of galaxy groups onto a massive z ~ 1 clus- 
ter accounts for a significant proportion (> 30%) of the galaxy 
population, and moreover, galaxies that derive from group en- 
vironments are more massive. Indeed, the two outliers are 
both more massive than the majority of main sequence galax- 
ies and exhibit suppressed SFRs. However, the sparse spec- 
troscopy at these high radii precludes us from reliably identi- 
fying infalling groups. 
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FIG. 9. — A caustic diagram for all spectroscopically confirmed cluster 
members. The symbols are the same as in Figure|7] We plot caustic profiles — 
lines of constant (r/r20o) X (Av/cr v ) — at 0. 1 and 0.4 to isolate virialized 
(within the inner caustic), backsplash (between caustics), and infall galaxies 
(along and outside the outer caustic). 



5.3. A Kinematic Approach to the Environmental 
Dependence on the Specific Star Formation Rate 

In hopes of disentangling radial contamination from envi- 
ronment, we consider the SSFR as a function of r x v, nor- 
malized by the cluster velocity dispersion and r 2 nn, in Fig- 
ure [TOl As shown in Figures|8]|9] and fig. 3 from lHaines et al.l 
( 120121) . low values of r x v, < 0.1(r 2 oo Xff„), should represent 
a more dynamically-isolated cluster core, as it is dominated 
by galaxies accreted at early times, which have spent more 
time in dense cluster regions. Similarly, galaxies recently ac- 
creted possess higher caustic values of > 0.4(r2oo x °~v), as 
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they are more likely to have high projected radii and veloc- 
ities. In Figure [TO] there is a clear segregation between the 
two star forming populations when using a caustic approach 
to define environment. The sub-main sequence galaxies (pink 
squares) are preferentially at low values of r x v and well con- 
tained within 0.4 (V200 x Ot>) with the exception of one galaxy; 
the main-sequence branch (cyan squares), however, does not 
penetrate below a value of 0.1 (r2oo X & v ), but rather is dis- 
persed throughout intermediate and high values. The isolation 
of the galaxies on the sub-main sequence yields a prominent 
depression in SSFR towards the "caustic core" of the clus- 
ter, namely, with galaxies that exhibit both low projected radii 
and low line-of-sight velocities. Splitting the galaxies into the 
aforementioned r x v bins (shown by vertical lines in Fig- 
ure [Tot and calculating the average SSFR of the MIPS cluster 
members in each bin (green stars) reveals a 0.9 dex decline be- 
tween the highest and lowest caustic bins. With this definition 
for environment, the cluster core is more of a signpost for the 
dynamic history of its constituents, as the time-averaged den- 
sity of a galaxy should increase with its time since accretion. 
If this dynamic core more accurately represents the virialized 
cluster core than clustercentric radius or density probes, the 
interpretation that the SSFR of star-forming galaxies lacks a 
dependence with environment warrants caution. Moreover, a 
rapid quenching timescale would no longer be a necessary re- 
quirement; if younger star-forming galaxies on the sub-main 
sequence are backsplash galaxies, we can use the dynamical 
timescale to place a lower limit of ~1.3 Gyr on the duration 
of quenching. 
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FIG. 10. — The SSFR for star-forming galaxies as a function of caustic 
environment, (r/r20o) X (Au/uu ). The sub-main sequence galaxies (pink 
squares) primarily inhabit low values of (r/r2on) X (Av/a v ) compare to 
the main-sequence galaxies. In contrast to Figure [4] there is now a clear 
depression (0.9 dex) of the average SSFR (green stars) towards low caustic 
values, where we expect to isolate galaxies accreted at earlier times. The gray 
vertical lines indicate the limits of the bins. 



MIPS cluster members form a young envelope of the 
cluster population; they span the same range of stel- 
lar masses and D4000 values as coeval field galaxies, 
suggesting they primarily belong to a recently accreted 
population. 

There is a bimodal or flat distribution of velocities for 
the MIPS members, peaking at ± 1000 km s -1 , which 
is roughly the cluster velocity dispersion. This sharply 
deviates from the older cluster members, which form 
the expected Gaussian distribution of velocities. There 
is only a 6% probability that these two distributions 
derive from the same parent population; therefore, the 
MIPS sources seem to represent a recently accreted, un- 
virialized population. 

When local environment is defined as projected clus- 
tercentric radius or density, it appears to have no effect 
on the average SSFR of star-forming galaxies (but see 
conclusions 5-8 for an alternative explanation). How- 
ever, when accounting for the total mass from all clus- 
ter members, there is a sharp decline in the total SSFR 
towards low radii, indicating that the SFR-density re- 
lation is already established in the densest regions at 
z ~ 0.9. 

There exist two branches of SFR/SSFR as a function 
of stellar mass — one that is consistent with the star- 
forming main sequence of increasing SFR with mass, 
while the other displays a flat trend with suppressed 
SFRs compared to the z ~ 0.9 field. This double- 
sequenced distribution appears to be unique to the clus- 
ter. Stacking the spectra of galaxies along each branch 
separately supports this by revealing distinct star forma- 
tion histories: the former are likely obscured galaxies 
experiencing either continuous or bursty star formation, 
while the latter lacks Balmer absorption, has a strong 
4000-A break, and weak [Oil] emission — features that 
are more typical of passive or quiescent galaxies. 



Utilizing the results from lHaines et al.l (120121) . we show 
that recently accreted galaxies contaminate all pro- 
jected radii, while the earliest accreted galaxies are pri- 
marily at low radii, with only a little spillover at larger 
radii. Assuming star forming galaxies are preferentially 
an infalling or recently accreted population, the star 
forming population is therefore inherently more likely 
to be a contaminant in projected environment, which 
could contribute to the flat trend between the SSFR of 
star-forming galaxies with environment. Therefore, low 
radial bins are contaminated by recently accreted MIPS 
sources that have not yet (necessarily) passed through 
the densest regions of the cluster or core and could aug- 
ment the SSFR. 



6. CONCLUSIONS 

We have presented an infrared study of SpARCS 
J16 13 14+564930, a z — 0.872 massive galaxy cluster com- 
prising 85 mass-limited, spectroscopically confirmed mem- 
bers, 16 of which have been detected at 24 /im. After com- 
bining the stellar ages, SFRs/SSFRs, and kinematics of dusty 
star forming galaxies, a clean narrative of the history of these 
cluster members begins to unfold. We summarize the main 
properties and trends of the MIPS galaxies as follows: 



A caustic digram of line-of-sight velocity versus pro- 
jected radius reveals that there is a mix of galaxies from 
each of the two star-forming sequences at any given ra- 
dius. We are able to successfully isolate galaxies that 
were accreted at early times from recently accreted ob- 
jects via caustic values of (r/r2oo) x (Av/cr„) < 0.1 
and > 0.4, respectively. Intermediate values contain a 
mix of accretion histories, and should also contain the 
statistical majority of backsplash galaxies. 
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7. Applying caustic profiles of constant (r/r2oo) x 
(Av/a v ) to our sample allows us to kinematically clas- 
sify MIPS galaxies as virialized, backsplash, or in- 
falling based on their age and SFR. Galaxies with lower 
star formation rates (i.e., on the sub-main sequence) 
appear to belong to either the virialized or backsplash 
populations. Further differentiation between these pop- 
ulations arises from stellar age, as measured through the 
depth of the 4000-A break. Although not definitively, 
older galaxies with D4000 > 1 .4, are likely to be viri- 
alized within the 0.1 (Y200 x &v) caustic profile. Back- 
splash galaxies appear to lie in intermediate regions of 
0.1 < r/7'200 x v/a v < 0.4. Infalling galaxies pri- 
marily occupy the region outside the 0.4 (Y200 x a v ) 
caustic profile and have SFRs consistent with the main 
sequence. 

8. Using a caustic definition for environment reveals al- 
most an order of magnitude (0.9 dex) depression of 
SSFR at low (r/r2oo) X (Av/<r v ) values, in contrast 
to the flat correlation observed with projected radius or 
density. This suggests that environment may suppress 
star formation as galaxies fall deeper into the cluster 
potential, with a minimum quenching timescale given 
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